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Abstract 

Consider the constitutive law for an isotropic elastic solid with the 
strain-energy function expanded up to the fourth order in the strain, 
and the stress up to the third order in the strain. The stress-strain re- 
lation can then be inverted to give the strain in terms of the stress with 
a view to considering the incompressible limit. For this purpose, use 
of the logarithmic strain tensor is of particular value. It enables the 
limiting values of all nine fourth-order elastic constants in the incom- 
pressible limit to be evaluated precisely and rigorously. In particular, 
it is explained why the three constants of fourth-order incompressible 
elasticity fi, A, and D are of the same order of magnitude. Several ex- 
amples of application of the results follow, including determination of 
the acoustoelastic coefficients in incompressible solids and the limiting 
values of the coefficients of nonlinearity for elastic wave propagation. 
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1 INTRODUCTION 



Extracting the condition for incompressibility from a stress-strain relation 
can be an ambiguous process because it leads to an infinite limit for one (or 
more) of the elastic stiffnesses and eventually, to the appearance of a hy- 
drostatic stress term proportional to an arbitrary Lagrange multiplier (to be 
determined from boundary and/or initial conditions) [H |2] . In linear isotropic 
elasticity, the relation between the infinitesimal stress cr and the infinitesimal 
strain e reads 

a = Atr(e)<5 + 2 fie, (1) 

where A and fi are the Lame coefficients and d is the identity. The incom- 
pressible limit is equivalent to the condition tr e = 0, which leads to the 
limiting case 

A — > oo, cr = — pd + 2 fie, (2) 

where p is an arbitrary scalar. By way of contrast, a strain-stress relation 
is more amenable to the imposition of incompressibility, because it leads to 
unambiguous, finite limit(s) for one (or several) compliance (s) [3J. Hence, in 
linear isotropic elasticity, we go from 

e = --tr( C r) c 7 + __< T , (3) 

where v is Poisson's ratio and E is Young's modulus, to the limiting case 

113 
v-±~, e = tr(cr)<5 H cr, (4) 

2' 2E v ; 2E ' v ; 

and there is no arbitrary quantity. 

Turning now to nonlinear isotropic elasticity, we must first of all make 
a choice of the measures of stress and of strain. Physicists and acousticians 
seem to favor the pair consisting of the Green-Lagrange strain tensor e and 
the second Piola-Kirchhoff stress i, and they expand, in the so-called weakly 
nonlinear theory, the strain-energy density W in terms of three isotropic 
invariants of the strain. Hence, the Landau and Lifshitz[3] expansion can be 
conducted to fourth-order as[5j 

W = ^II + fih + p 3 + Bhh + jlf + Ehh + Fllh + Gil + Hlf, (5) 
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where A, B, C are the three third-order constants, E, E, G, H are the four 
fourth-order constants, and 

Ji=tr(e), / 2 =tr(e 2 ), 7 3 = tr(e 3 ), (6) 

are respectively of first, second, and third order in the strain. (Note that we 
have placed overbars on the elastic constants so as to avoid conflict with the 
standard notation E for Young's modulus, which is used extensively in what 
follows.) Thus, terms in l\ and I2 are "second-order" terms; terms in If, 
I1I2, I3 are "third-order" terms, and so on. However, the incompressibility 
limit is not easily implemented with this choice, because it is a combination 
of the invariants [6] , specifically 

h - h + Ti + m - 2/i/ 2 + 1/ 3 = 0. (7) 

This, in particular, means that I± is now a second-order quantity. This makes 
it complicated to arrive at the fourth-order expansion of incompressible non- 
linear elasticity, for which the strain-energy function has the form 

W = fil 2 + p 3 + Dll (8) 

where D is another constant. 

However, it is advantageous to use, instead of Green-Lagrange strain, 
an alternative measure of strain, namely the logarithmic strain, since this 
leads to a very simple expression of the incompressibility constraint. The 
logarithmic strain e is defined by 

e=iln(5 + 2e), (9) 

and we consider the three independent invariants 

h = tr{e), J 2 = tr(e 2 ), J 3 = tr(e 3 ), (10) 

respectively of orders one, two, and three, analogously to Ji, J 2 , I3. Then 
incompressibility is expressed exactly as 

h = 0, (11) 
which must hold identically for all stresses and strains. 
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These considerations suggest the following protocol for finding the lim- 
iting values of the elastic constants in incompressible weakly nonlinear elas- 
ticity (Sections |2] and |3]). First, write down the relation between stress and 
strain. Since we are considering isotropic materials, the relevant stress mea- 
sure t conjugate to logarithmic strain e is the Kirchhoff stress tensor on 
rotated coordinates, related to i by t = (8 + 2e)i (for a general discussion of 
conjugate stress and strain tensors see OgdenJT]). Then invert this relation 
to find the strain in terms of the stress, and, more particularly, I\ in terms of 
invariants of t. Finally, note that l\ must be zero for all t. In what follows 
we conduct this process explicitly for third-order elasticity. Then we present 
the results of the fourth-order case, omitting many of the (cumbersome and 
lengthy) details of the calculations. In Section HJ we give a few examples 
of applications, and we see that results established in compressible elasticity 
can be taken to their incompressible limit, without having to re- write and to 
re-solve the equations of motion and the boundary conditions. 

Previously, Hamilton et al. [8] have shown that there should remain only 
three elastic constants in the incompressible limit of Eq. §5§ (OgdenjH] had 
in fact proved this result 30 years earlier). They also found some limited 
information on the behavior of the other six constants, by extending the 
work of Kostek et al. [10] from third-order to fourth-order elasticity. Based 
on a comparison with the equation of state of inviscid fluids, it was found 
that /i and A remain, a new fourth-order constant D emerges, and 

A -»■ oo, B = -A, E 

F = -C, G = ^, D 

see Zabolotskaya et al. [TT] and Jacob et al. [5] for the latter identity. However, 
the behavior of C, F, and H remains undetermined. Moreover, the full 
comparison with fluids ultimately leads to the identities: fi = 0, A = 0, 
D — 0, which, in considering an elastic solid, are not satisfactory because 
they imply that some information might be missing from the list ( IT21) . where 
/i, A, and D could play a role. One can only go so far in comparing the 
behavior of a solid to that of a fluid: a particular disconnect emerges when 
comparing the behavior of a solid in the incompressible limit, where the speed 
of a longitudinal wave should tend to infinity, to that of an isentropic liquid, 
where the speed of a longitudinal wave is finite (Section I4.4p . The differences 
existing between a compressible solid and its incompressible counterpart must 
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be tackled within the framework of solid mechanics. 

In effect, many questions have remained open in the literature, as attested 
by the following comments. Domariski[T2"] remarks that "the details of the 
derivation [of Eq. (jSJ)] are not quite clear from the mathematical point of 
view", and that "surprisingly, the experiments confirm that, in spite of being 
a combination of the higher order constants, the fourth-order constant D 
is of a similar order of magnitude as the second-order shear Lame constant 
/i and the third-order Landau constant A" . Catheline et al. [13] measure A 
and B for Agar-gelatin based phantoms and note that "the huge difference 
between these third-order moduli is striking since in more conventional media 
such as metal, rocks, or crystals they are of the same order." They then 
propose an "intuitive justification" for this difference, which "does not hold 
for a theoretical explanation" . Jacob et al. [6] record that "surprisingly, the 
experiments confirm that, although expressed as function of compression 
moduli A, B, G, the moduli D is of the order of magnitude of the so-called 
shear moduli ft and A" and add that "no explanation has been given for this 
order of magnitude." We address each of these points in the course of this 
paper. 



2 THIRD-ORDER INCOMPRESSIBILITY 

In this section we focus on third-order elasticity, so that the strain-energy 
function (jSJ) reduces to 

W = ^II + fih + p 3 + BI x h + ji?- (13) 
At this order, the inversion of Eq. (j9j) gives e = e + e 2 + 2e 3 /3, so that 

h = h + h + \h, / 2 = /2 + 2/ 3 , h = h- (14) 

It follows that the third-order expansion of W in terms of the invariants of 
the logarithmic strain e reads 

W = + nh + h + (B + \)IJ 2 + jlf, (15) 

which is also written as 

W = — ^ -I? + , E A 2 + — 1 3 + Bhh + -If, (16) 

2(1 + i/)(1-2i/) 1 2(l + i/) 3 1 2 3 J ' y J 
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where E, v are second-order constants (Young's modulus and Poisson's ratio) 
and A, B, C are third-order constants with respect to the logarithmic strain, 
with the connections 



^ 3A + 2/i A 

E = — U, V = 77 T, 

A + /x 2(A + /x)' 

A=A + Q^ B = B + X, C = C 



(17) 



to the Lame and Landau coefficients. The (conjugate) stress t = dW/de 
then expands as 



Eu 



(1 + i/)(1-2i/) 



h + Bh + Clf S 



E 



l + v 



+ 2Bh e + Ae\ (11 



Now we introduce the stress invariants T\ = tr(t) and T 2 = tr(t 2 ), and it 
follows by taking in turn the trace of Eq. ( TT8|) and then of its square that 



To 



T 2 



E 



1 - 2v 
E 2 

{l + v) 2 

E 2 
(l-2u) 



h + (2B + 3C)lf + (A + 3B)I 2 , 



"(2 -v)j2, r 



1 - 2i/) 2 

2 



(19) 



correct to the second order in the strain. These can be inverted, to give If, 
I2, and then 

/i = «T 1 + /3T 2 + 7 T 1 2 , (20) 
where the constants a, (3, 7 are given by 



a 



\-2v 
~~E~' 



(3 = -(A + 3B) 



:i + ^r(i-2^) 

E 3 



^(^3 B ) ^-l 1 - 2 -» -(2 B + 3 C )^. (2!) 



£ 3 



£ 3 



For the exact incompressibility condition Eq. (fTTj) to hold for all stresses 
and strains, we must have a = = 7 = 0. First, a = is clearly equivalent 
to either of 

1/ — ► 1/2, A -)■ 00, (22) 
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whilst both E and \x remain finite, with [i — > E/3. 

Second, note that the stress must remain finite in the incompressible limit. 
Hence, the last term in the expression for t in Eq. f fT8|) remains finite, so 
that A remains finite: 

A/fi = 0(l). (23) 

With these two conditions for incompressibility, the expression for /3 in 
Eq. ( 12"T|) reduces to 

27g(l - 2u) 

P 4£ 3 

Notice that the first term in the expression for t in Eq. (jl8p must also remain 
finite, which means that (1 — 2v)~ 1 I\ remains finite. Clearly, the conditions 
"/3 = 0" and "(1 — 2u)^ 1 (3 is finite" are fulfilled simultaneously when B 
remains finite: 

B/n = 0(l). (24) 

With the three conditions for incompressibility Eqs. (l2"2"]) - (124j) . the ex- 
pression for 7 in Eq. fj2Tj) reduces to 

7=-?^. (25) 

Recall that (1 — 2z/)~ 1 i 1 and hence, (1 — 2z/) _1 7 remain finite. It follows that 

(l-2z/) 3 C^0, C/fi= 0{X 2 /fi 2 ). (26) 

In summary, the conditions given in Eqs. (122"]). fl23|) . ( pij) . and (j26"j) are 
necessary and sufficient for incompressibility. 

For an incompressible material (with I x = 0) the strain energy W there- 
fore reduces to 

W = jl 2 + p 3 , (27) 
or equivalently, in terms of the invariants of the Green strain tensor, 

W = |(/ 2 - 2/ 3 ) + jl 3 = fj 2 + p 3 , (28) 

where we have used Eqs. ffT4"]) - (fTT|) . We use Eqs. ffT7|) to find that, in terms of 
Lame and Landau constants, the incompressible limits are Eqs. (T2"2"|) together 
with 

A = A-6^i = fiO{l), B = B-\ = XO{l), C = C = fiO{X 2 /fi 2 ), 
(1 - 2u)B -E/3, (1 - 2ufC -> 0. (29) 
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Notice that C varies in a quadratic manner with respect to A as it goes to 
infinity in incompressible solids, not in a linear manner, as is incorrectly 
reported in Refs.JHl [HI [15]; see also Section I4.4[ where it is shown that if C 
were linear in A then longitudinal waves would propagate with finite speed 
in the incompressible limit. 

Bearing in mind the advantage of the strain-stress relation in the 
linear theory we now take note of its counterpart for the present theory. On 
use of the definitions of 7\ and T 2 and Eq. (fT9|) we may, after lengthy algebra, 
invert the stress-strain equation (fl8|) to give 

e = —^t-—T x 6 — 1 + — — — \vA - (1 - 2v)B\ {2T x t + T 2 5) 

+ \ZBv{2 - i/)(l - 2v) - 3Au 2 - C(l - 2z/) 3 l Tf*. (30) 
In the incompressible limit as embodied in Eqs. f l22|) - f[26l) . this reduces to 

e = ^ (2£ 2 + 3^) (t - |T 1( 5) - ^ (* 2 - |T 2 <5) , (31) 

from which it follows immediately that tr(e) = 0. 

Conversely, the stress-strain relation must accommodate the internal con- 
straint of incompressibility by introducing a Lagrange multiplier, denoted p, 
in the following expansions, 
2E 

t = -pS + — e + Ae 2 , t = -p(8 + 2ey 1 + 2{ie + Ae 2 , (32) 
3 

for the stresses, where we recall that for an isotropic material, t = (6 + 2e)i. 



3 FOURTH-ORDER INCOMPRESSIBILITY 

We now extend the above analysis to include fourth-order terms in the strain- 
energy function, as in (jSJ), and we work in terms of the logarithmic strain 
tensor and its invariants. We use the fourth-order expansion e = e + e 2 + 
2e 3 /3 + e 4 /3 and the identity tr(e 4 ) = (l/6)/ 1 4 -/ 1 2 / 2 + (l/2)/f+ (4/3)/iJ 3 
to establish the following connections between the invariants of e and those 
of e: 

h — h + h + \h + — 3-^1 ^2 + + 1^1 > 

I 2 = h + 2I ?J - \l\h + fhh + |if + M 4 , 

7 3 = h + 4/i/ 3 - 3/ 2 / 2 + |/ 2 2 + |/* (33) 
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The strain energy may then be expressed as a function of the invariants 
of e, in the form 

W = ±1* + iih + jh + Bhh + |i? + Bhh + Flih + Gil + ml (34) 
where, in addition to the connections (117]) . the elastic constants are given by 

„ - - 4A 28u 2A - - , 7/i 

£ = £ + 25 + — + — + — , J r = F + C-A-— , 

3 9 3' _ 3' 

Note that it is a simple matter to invert these relations using Eqs. (|17|) . to 
give 

The stress t, to third order in e, is then given by 

t = XhS + 2//e + B(I 2 5 + 2he) + Ae 2 + Ci?6 

+ £(I 3 d + 3he 2 ) + 2T{hI 2 6 + lie) + 4QI 2 e + mifd. (37) 

By defining T\ = tr(t), T2 = tr(i 2 ), T3 = tr(i 3 ), similarly to the previous 
section, by manipulating the equations to third order, and after some consid- 
erable algebra (which is omitted), we obtain an extension of the expansion 
(JH]) to give 

h = olT x + (3T 2 + 7T 2 + a% + 0%T 2 + 7'ri 5 , (38) 
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where a, fl, 7 are as given by (}2"Tj) . while a', (3', 7' are defined via 



E 5 a' 



(1 + z/) 3 (l - 2u) [2(1 + i/).A(.A + 3B) - 3E£] , 

(1 + u) 2 {l - 2u)^2{A + 3B) [(1 - 2i/)(5 - 4i/)S + 3(1 - 2v) 2 C 

-i/(4 + - £ [3(1 - 5z/)£ + 2(1 - 2i/)(3J 7 + 20)] }, 



E 5 (3' 



ES = 2v 2 {2-v){l-2v) 2 {A + 3B) 2 

-6i/(l - i/)(l - 2z/) 3 (^ + 3B){2B + 3C) 

+ 2(l-2z/) 5 (2£ + 3C) 2 

+3z/ 2 (l - 2i/) [2(1 + i/).A(.A + 3B) - 3E£] 

- i/(2 - z/)(l - 2z/) 2 [6(l + i/)(.A + 3£)£ - £(3£ + 6 7 + 40)] 

+ 2(1 - 2v)\\ +u)(A+ 3B)C - E{7 + QU)}. (39) 



The exact incompressibility constraint Eq. All]) is enforced for all stresses 
when a = /3 = 7 = and a' = = 7' = 0. The first three conditions lead 
to the limits Eqs. (121Z1) - ([2"6"1) . thereby bringing great simplifications in the 
expressions above. Hence in the incompressible limit we have 



which must tend to zero as v — > 1/2. However, on inspection of ( )37|) it can 
be seen immediately that for the stress to remain finite in the incompressible 
limit we must have that £ is finite: 



and a' above does indeed tend to zero as v — > 1/2. As in the previous section, 
we see that the first term in Eq. (1371) remains finite in the incompressible 
limit when (1 — 2v)~ x I\ remains finite. Clearly here, (1 — 2v)^ x ol remains 
finite. 

Equally, for the term in Q in Eq. fl37j) to remain finite, Q must itself be 
finite in the limit: 



E A oc 




2u)£, 



(40) 



e//i = o{i), 



(41) 




(42) 



E A (3' = -f {\-2vfT. 



(43) 
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For f} 1 to tend to zero, and (1 — 2v) 1 j3' to remain finite, we must enforce the 
following behavior for J 7 : 

(l-2z/)lF^0, F/yL = 0{\/yL). (44) 

In these limits, the /i/ 2 J-"-term in the expression Eq. f p?T|) for the stress 
remains finite, whilst the if .F-term vanishes. 

It remains to consider %. Using the limits above, we see that 7' behaves 

as 

ES = -12(1 - 2vfU. (45) 
For this to tend to zero, and (1 — 2v)~ lr y' to remain finite, we require 

(l-2z/) 4 ^^0, n/fi = 0(\ 3 /fi 3 ). (46) 



We may then check that the limiting value of the last term in Eq. f l37|) is 
finite. 

In summary, we must have for the fourth-order constants associated with 
the logarithmic strain 

e/n = o{i), ^/ai = o(a/a*), g/n = o{i), H/» = o(\ 3 / f i 3 ), 

(47) 

which are necessary and sufficient for incompressibility at this order. 

For the fourth-order constants associated with the Green strain, we use 
Eqs. (J36J to find 

E/n = 0(\/fj.), F//i=0(A 2 //i 2 ), G/fx = (D(\/fx), H/fi = 0(A 3 //i 3 ), 

(48) 

and, specifically, 

- AE E 
{ l_2u)E^— 1 (i-2v)G^-. (49) 

6 

For incompressible solids, 7i = for all stresses, and the fourth-order 
expansion Eq. (1341) reduces to 



w = fii 2 + ji 3 + gil (50) 

and only three constants remain[9]. Equivalently, in terms of the invariants 

of the Green strain tensor, the expansion of W reads as Eq. (JSJ), where D is 
defined by Eqs. ( TT2l) or, equivalently, by 

D = g-^ + ^ (si) 

2 6' V ; 
11 



making it explicit that it is of the same order as [i (and thus, as A): 

D/n = OQ). (52) 

Turning our attention to the stress, we see that all the terms multiplying 
8 in fl37|) are absorbed by the arbitrary hydrostatic stress, to give, in the 
limit, 

t = - P S + 2/ie + Ae 2 + AGI 2 e. (53) 

By the Cayley-Hamilton theorem we have (for an incompressible material, 
where I\ — 0) 

e 3 = |/ 2 e + |/ 3 <5, (54) 

and hence we have 

t = -pd + 2fie + Ae 2 + Ag(2e 3 -p 3 S), (55) 

and finally, by adjusting the hydrostatic term by introducing p' = p+8QI 3 /3, 

t = -p'S + 2fie + Ae 2 + 8£e 3 , (56) 

where only the three constants of fourth-order incompressible elasticity ap- 
pear, and no invariant. We recall that the corresponding measure of stress 
conjugate to the Green strain is given by the connection i = (6 + 2e)~ 1 t, 
which applies for an isotropic material, yielding 

t = -p\S + 2e)- 1 + 2/ie + Ae 2 + 8De 3 . (57) 

Note that the Lagrange multiplier p' must figure in the expressions for the 
stress (see also Eq. (|32|) ). but has been omitted in the expression for the 
stress in several paperspTJ EJ [16l IT?] . 

The counterpart of the strain-stress relation ( 1301) for the fourth order is 
very lengthy and is not written here. We note, however, that the contribution 
to e additional to the first and second-order terms in fl30|) for the third order 
in the stress has the structure 

ai t 3 + a 2 T 1 t 2 + (a 3 T 2 + a 4 T 2 )t + (a 5 T 3 + clqT^ + a 7 T?)5, (58) 

where a 1; a 2 , . . . , a 7 are constants that are collectively functions of /i, A (or 
E, u), A, B, . . . , %. With the exception of a± = A 2 /16fi 5 these expressions 
are very lengthy and therefore omitted. However, in the incompressible limit 
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the strain-stress relation simplifies substantially and the extension of (1311) 
becomes 



e = 




^-{2E 2 + ?,AT 1 ){t- 



\TiS) 



27 



(T? -3T 2 )(8EQ -A 2 )(t 



ins). (59) 



32E 5 



As in the third-order case it is seen immediately that tr(e) = 0. In deriving 
(|59p we have used the Cayley-Hamilton theorem for t to eliminate t 3 . 



4.1 Acoustoelasticity of bulk acoustic waves 

Hughes and Kelly [IS] used the acousto-elastic effect to evaluate experimen- 
tally the third-order elasticity constants, by measuring the speed of infinites- 
imal bulk homogeneous plane waves propagating in a solid subject to a small 
pre-stress. We summarize their results in Table I, where we use the layout 
of Norris[l9] and the definition 



of the acousto-elastic coefficient K, where p is the mass density in the un- 
strained state, and a is the pre-stress (which is either hydrostatic or uniaxial). 

The incompressible counterparts to these formulas are readily established 
by use of Eqs. f )22|) and fl29|) . They appear in the last column of the table. 
In particular, it is seen that the speed of longitudinal waves is infinite, as 
expected in an incompressible solid. The formulas for shear waves are in 
accord with those established using different means by Gennisson et al. |20] 
and by Destrade et al. [21] in the case of uniaxial pre-stress. For hydrostatic 
pre-stress, we note that the speeds of the shear waves are unaffected by 
the hydrostatic stress, which is also to be expected for an incompressible 
material. 

4.2 Acoustoelasticity of surface acoustic waves 

Hayes and Rivlin [22] computed the acousto-elastic coefficient for surface 
acoustic wave (SAW) propagation. For a wave propagating in the direction 



4 EXAMPLES 
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Table I. Acousto- elastic coefficients for bulk acoustic waves in compressible and incompressible solids (here 

3k = 3A + 2/ij. 
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_L stress 


_L stress 
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of uniaxial pre-stress of (small) magnitude a, it is defined by Eq. (ISO]) where 
v is now the SAW speed. Using the results of Tanuma and Man[23, 24J, we 
present it in the form 

K = 1 - a7 22 - &723 - C733 - G?744 (61) 

for a compressible isotropic elastic solid, where the 7's are defined in terms 
of the Lame constants by 

722 = (A + 2/i) [-8(A + y) + 2(5A + 6/i)X - (2A + 3/i)X 2 ] /A, 

723 = 4A(1 - X) [4(A + /i) - (A + 2/i)X] /A, 
733 = [1 - /iX/(A + 2/i)] /A, 

744 = -8(A + 2/i - 11X) [2(A + n) - (A + 2//)X] /A, 

A = (A + //) [8(3A + 4//) - 16(A + 2/i)X + 3(A + 2/i)X 2 ] , (62) 

the non-dimensional quantity X is the unique real positive root to Rayleigh's 
cubic|25 



X 3 - 8X 2 + 8 3 4±^X - 16±±JL = 0, (63) 
A + 2/i A + 2/i v ; 

and a, b, c, d depend on the second- and third-order constants according to 
_ _ (4A + 3/i)(A + 2/i) + 2(A + p)A 2(2A + 3p)B ^ 2C 



(3A + 2/i)/i (3A + 2/i)/i (3A + 2/i)/i 3A + 2/i 

(A + /i)A (A + 2fi)B 2C 



(3A + 2/i)/i (3A + 2/i)/i 3A + 2/i ' 
(2A + /i)(A + 2/i) \A 2(X-fj,)B 2C 



(3A + 2/i)/i (3A + 2/i)/i (3A + 2/i)/i 3A + 2/i' 

d= A + ^ + ( A + 2 ^ + g (64) 
(3A + 2/i)/i 4(3A + 2/i)/i 3A + 2/i' v ; 

Substantial simplifications occur when the incompressible limits (I22 p - (|29[) 
apply. In particular, X is now a definite number [25], namely X = 0.9126, 
the root of X 3 - 8X 2 + 24X - 16 = 0, and 722 = -723/2 = 733 = 2(X - 
1)(4 — X) / (26 — 16X + 3X 2 ). Because of these latter relationships, the limit 
of K remains finite since, even though each of a, b, and c goes to infinity 
as 0(A//i), the combination a — 2b + c tends to a finite limit, specifically 
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a — 2b + c — > A/(3p), whilst d — > A/(12p). The final result is that the 
acousto- elastic coefficient for SAWs in incompressible media is 

A 

K = 1 + 0.9126 , (65) 

as established differently by Destrade et al.\21\. We emphasize that this result 
applies for the situation in which the stress is uniaxial. It is interesting to 
note in passing that in the corresponding plane strain problem (in the (1,2) 
plane with stress a in the Xi direction) it can be shown that K = 1 — X/2, 
which is independent of the third-order constant A. 



4.3 Solitary waves in rods 



Porubov[26j showed that the propagation of nonlinear strain waves in an elas- 
tic rod with free lateral surface is governed by the so-called double- dispersive 
equation (DDE). For solids with the third-order strain energy density Eq. 
([13]), the DDE is 



v tt - aiv xx - a 2 



(v 2 ) xx - a 3 v xxtt + a 4 v xxxx = 0, (66) 



for a strain function v = v(x, t), where x is the space variable in the direction 
of propagation, t is time, subscripts denote partial differentiation, and 

E p vjy-1) p2 uER 2 

on = — , "2 = ct 3 = 7, R > "4 = 7, — • (67) 

p 2p 2 2p 

Here E is Young's modulus, v is Poisson's ratio, R is the rod radius, p is the 
mass density, and (3 is the nonlinear parameter. (Note that this is different 
from the defined in Eqs. (I2T]) .) Explicitly, 

P = 3E + (1- 2vf(B + C) + 2(1 - 2z/)(l + u)(A + 2B) + Qu 2 A. (68) 

Now, in the incompressible limit described by Eqs. (!22|) - (!29l) . this non- 
linear parameter reduces to the very simple expression 

/3 = |(4/i + A), (69) 

and the a's simplify accordingly. The analysis of Porubov can then be carried 
through, in particular to study solitary waves and solitons. It should be 
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pointed out that the sign of (3 can be determined for all incompressible third- 
order solids. Indeed, it is well-known that, at the same level of approximation, 
the model described by the Mooney-Rivlin strain energy density 

W = C 10 [tr(C) - 3] + CbilMC- 1 ) - 3], (70) 

where C = 2e + d is the right Cauchy-Green strain tensor and C\o and Cqi 
are constants, is equivalent to the strain energy density Eq. fl28|) . with the 
connections 

fi = 2(C 10 + C 01 ), A = -8(C 10 + 2C 01 ). (71) 

Clearly, it follows that here we have (3 = — 32C*oi. However, it is also well- 
known that the governing equations of motion for Mooney-Rivlin solids are 
strongly elliptic when [27] Cqi > 0. Provided that this condition is satisfied, 
we deduce that 

< 0, (72) 

for all incompressible third-order solids and therefore that a compressive 
solitary wave emerges from an initial localized compressive input [26J. 

In fact, Porubov and Maugin[28] show that fourth-order elasticity is re- 
quired to account for the possibility of simultaneous compressive and tensile 
solitary waves. They use the Murnaghan|29j counterpart to Eq. ([5]), where 
the expansion is carried out in terms of the principal invariants of e, which 
we write as 

ii = tr e, i 2 = \ [(tr e) 2 - tr(e 2 )] , i 3 = det e. (73) 
The Murnaghan strain-energy function is then given by 

A + 2u . 2 I + 2m .o . . . , 4 .o. . . .o 
W = — - — i l - 2/ii 2 H — 1\ - 2mi x i 2 + ni 3 + V\% x + v 2 i x %2 + v 3 iii 3 + z/ 4 z 2 , 

(74) 

where m, z/ 4 are constants. The correspondence between the Landau 

and the Murnaghan constants is easy to establish. We find the connections 

A 

m = — + B, n = A, l = B + C, 

Ul = E + F + G + H, v 2 = -3E -2F- AG, u 3 = 3E, z/ 4 = AG, 

(75) 
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and thus, the incompressible limits 



m/fi = 0{X/fi), m/fi = 0(1), l/fi = 0(X 2 /fi 2 ), 
v 1 /n = 0(\ 3 /n 3 ), u 2 /fi = 0(X 2 /fi 2 ), 

v 3 /fi = 0(X/fx), vt/n = 0{\/fi), (76) 

with the specific limits 

(1 - 2u)m -> -E/3, (1 - 2v)v 3 -> 4E, (1 - 2i/)i/ 4 ->■ 2£/3. (77) 

Now, when fourth-order terms are taken into account, the corresponding 
DDE has an extra term[28| and becomes 



v t t ~ u\v xx - « 2 {v 2 ) xx - a 3 v xxtt + a 4 v xxxx - a 5 (v 3 ) xx = 0, (78) 

where ati, a 2 , a 3 , a 4 are still given by Eqs. (JB"7j) . and a 5 = 7/(3p), with 7 
(different from the 7 in Section [2J given by 

E-f = E 2 - 8Z 2 (1 - 2z/) 5 (l + 1/) - 32m V(l - 2u)(l + z/) 3 
- 8nV(l - 2i/)(l + 1/) + 4/(1 - 2z/) 3 {£ - 4i/(l + i/)[2m(l + 1/) - ra]} 
+ 8m(£ + 4nz/ 2 )(l - 2i/)(l + z/) 2 + 12nEv 2 + 8v x E(l - 2vf 
- 8u 2 E{\ - 2u) 2 {2 - v)v + 8v 3 E(l - 2u)u 2 + 8z/ 4 £(2 - vfv 2 . (79) 

Clearly, the limits Eqs. (1761) — (1771) do not give a definite incompressible limit 
for 7. In particular, the terms proportional to (1 — 2u)m 2 and to ^4 tend 
to infinity, and the limits of the terms proportional to (1 — 2v) 3 lm and to 
(1 — 2v) 2 v 2 are not known. When 7 is written in terms of the Landau 
constants (using the inverse of Eq. (175]) ). similar ambiguities arise. However, 
if the expression for 7 is further transformed in terms of A, B, . . . , *H, the 
nonlinear constants associated with the logarithmic strain, then it is a simple 
matter to find its unequivocal incompressible limit. In the end we find that 
7 tends to 

7 = — - 3A + 180, (80) 

or equivalently, in terms of the incompressible Landau constant A and the 
constant D defined in ([12]) . 

7 = 14/1 + 6^ + 180. (81) 
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4.4 Nonlinear plane waves 



Wochner et al. [16j [30] also use the notation (3 to identify the coefficient of 
cubic nonlinearity for shear waves. Again, this is different from that used 
previously in this paper, and is given by 



A + 2/i + A + 2B + 2G - 



(A + 2/X + A/2 + Bj' 
A + /i 



(82) 



According to the first, second, and fifth limits in Eqs. ( TT2|) . this quantity 
should collapse to ft = 3(2/i + A)/(4/i) for incompressible solids. In fact, the 
true limit is 116 



ft 



1 + 



A + 2D 
2^i 



(83) 



because D, as defined in Eq. (fT2"j) . should remain finite, although until now, 
this behavior had not been proved rigorously. An alternative means of finding 
the correct limit is to rewrite the expression for ft in terms of the constants 
associated with the logarithmic strain, as 







3(l + i/) 
2E 



2G E 2(1 + !/)(! -21/) (A E 
y 6(l + i/) E V 2 2(l + z/) 



(84) 

Then the limit is clearly (3 = 9(2Q — E/9)/(4E), which is the same as Eq. 

Note that Wochner et al. [16] also provide fti, the Gol'dberg[3T] coefficient 
of nonlinearity for longitudinal waves, as 



A + 3B + C 
A + 2/i 



(85) 



If, as is incorrectly reported in Refs. [HI [HI [15] , C were to behave as C j \i = 
0(X/fj,), then (3i would remain finite in the incompressible limit, suggesting 
that longitudinal homogeneous plane waves are possible. As we have estab- 
lished in fl2"9"j) . C j ii behaves as C(A 2 //x 2 ), and ft\ therefore blows up in the 
incompressible limit, as should be expected, thereby precluding the existence 
of such waves. 



19 



5 CONCLUDING REMARKS 



Using the logarithmic strain measure, we are able to determine the exact 
behavior of the elastic constants of second, third, and fourth orders in the 
incompressible limit, as collected below. For the second-order Lame constants 

A -> oo, fi ->■ E/3, (86) 

as is well known; for third-order Landau constants, 

A/fi = 0(l), B/fi = 0(\/n), C//i = (9(A 2 //i 2 ), (87) 

and for the fourth-order Landau constants, 

E/vl = 0(\/h), F//i = 0(A 2 //i 2 ), 

G//i = 0(A//i), tf//i = 0(A 3 //i 3 ). (88) 

For the constants which vary linearly with A as it goes to infinity, the specific 
limits are 

(1 - 2u)B -»> -E/3, (1 - 2u)E -»■ AE/3, (1 - 2v)G -»• E/6, (89) 

as the Poisson ratio v — > 1/2. 

We have used these limits to show that it is easy to take known results 
of elastic wave propagation in compressible materials to the corresponding 
incompressible limits. In fact, other types of internal constraints could be 
accounted for just as easily[32, 33J. 

We conclude the paper with two remarks. The first is technical: it should 
not be forgotten that the hydrostatic term in the stress-logarithmic strain 
relation of an incompressible solid is an arbitrary Lagrange multiplier, to be 
determined from initial/boundary conditions. This Lagrange multiplier also 
appears in the stress-Green strain relation, see Eqs. (132]) and fl57|) . Omission 
of this term would therefore lead to incorrect solutions of the equations. 
The second is semantic: one gets the impression from reading the acoustics 
literature that 'soft' and 'incompressible' are two interchangeable adjectives. 
It should be clear that they are not, and Nature and Engineering provide 
many examples of hard materials which are incompressible (such as fully 
saturated soils in undrained conditions [31]) and of soft materials which are 
compressible (such as polyurethane foams [35J). 
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